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Abstract 

We discuss the two photon coupling of the lightest scalar meson on the basis of an extension 
of yPT. Using low-energy data on the pion form-factor and the 77 —>• 7r + 7r _ (7r°7r°) cross- 
sections as inputs, we find Y)a —> 77) = 0.126 keV. The smallness of the result and the 
relative weight between its components, < 1, suggests that the scalar 0 ++ meson 

is mainly a Q 2 Q 2 state. 

PACS numbers: 13.20.Jf, 11.30.Rd, 12.39.Fe, 12.38.Aw 


Motivation: The plethora of scalar mesons in QCD has a long and puzzling history. Probably, 
due to its elusiveness, the most interesting state is the isoscalar I = 0, er(600) [Tj. It is well known 
as a broad enhancement in very low-energy s-wave meson-meson scattering. The quark or gluon 
content of the cr(600) is not fully understood and the proliferation of models with seemingly 
different conclusions is disturbing m [ 3 ]. At the same time its underlying structure is a corner 
stone in understanding the realization of the mechanism for chiral symmetry breaking. 

In this paper we find indications of the Q 2 Q 2 content for the a state and estimate it by 
studying the processes 77 —> nir and the pion vector form-factor. The tetraquark structure of 
the lightest scalar was proposed long time ago owing to a possible strong diquark correlation g|. 
Our working framework runs in parallel to that in [5] with the only difference that we interpret 
their Lagrangian in an effective perspective by providing a counting power to the singlet field 0. 
Our main result is based on the comparison of two terms: the first one, already studied in 7], 
is given by the rescattering effects contribution to the 77 —>• 7T7 t —>■ Si decay and the second by 
the direct 77 -7 S 1 coupling. At the fundamental level the two-photon coupling for a generic S 1 
scalar meson is given by 

L = -|^c l7 Si. (1) 

There are many ways to couple the scalar singlet to the vacuum. If one considers that the 
spontaneously breaking of scale invariance is mediated via the trace of the energy-momentum 
tensor the coupling ci 7 is related to the scalar decay constant via the relation 



a <r(e + e —> hadrons) 
67 t <r(e+e~ n + ^~) 


and (0|6U|Si) 


—Mj^Fs, 


( 2 ) 


being 9 ^ the trace of the energy momentum tensor 

:= e g + 9 q = \p{a s )Gl v G» va + £ rmi>i (1 + 7 m(a B ))^ . 
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(3) 


Instead, if we consider that the scalar meson is a S-wave bound state of diquark-antidiquark pair 
the corresponding interpolating field can be constructed as 


jcr €abc€dec(n a C'ysd^i^Ud'YsGd e ), 


(4) 
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where latin indices denote color and C stands for the charge conjugation matrix. In the above 
expression the diquark is taken to be a spin zero color antitriplet and flavor antitriplet [5]. Then 
the coupling to the vacuum is given by m 

{0\3a\Si) = ~V2M 4 Si F Sl . ( 5 ) 


Setting the scheme: Let us first recall the main ingredients of the theoretical set-up. We 
shall consider an effective approach to QCD with two flavors in the isospin limit. The smallness 
of the values of the light-quark masses and the external momenta set a perturbative scheme out 
of the chiral symmetry limit. We count the pion and scalar field as O(p 0 ), derivatives, vector and 
axial-vector external currents as 0(p) and the scalar, pseudo-scalar external currents and scalar 
mass as 0(p 2 ). With this counting the leading order Lagrangian reduces to that presented in [TT 

£ 2 [0 ++ ] = + Fc ld S r + c 2d S 2 + • • •) <ut «**> 

+ + FClmSl + C 2rn Sf H-^ ((X+) — <X f + X)) ■ (6) 

The pseudo-scalar field is parametrized by the unitary matrix u{x) 2 = U(x ) = . 

Here F is the pion decay constant (F ~ 93 MeV), the fa's are fields for the pseudo-scalar Goldstone 
mesons and <Ji are the Pauli matrices. The basic building blocks are defined as 

X = 2 B 0 (s + ip ), 

= iv)DfJJu* = —iuDfJJ^u = u t , 

X+ = M t XM t + «X f w • (?) 

Next-to-leading order corrections, 0(p 4 ), come either through one-loo p g raphs or by higher 
order operators. In particular, the terms relevant to our study are explicitljO 

6 2 
A[0 ++ ] = Y, + ZiMfa'U + X U ] ) + Z 2 M 2 (D^UD>*U') - ^c^S.F^F^ , ( 8 ) 

i =5 

where M a stands for the singlet mass in the chiral limit and 

p 5 = , p 6 = \{ft v K,u v ]) , (9) 

with = uF R v v) ± v)F^u . The field strength tensors F^ R are related to the non-abelian 
external fields m- One salient feature of the field theory approach presented above is that it 
allows to separate between the direct and rescattering 77 couplings in a crystal clear fashion. The 
reason being that the direct process always involves the gauge invariant operator in ©, SiF^F^. 
Such separation is not always feasible using dispersion relations. 

Charged pion-pair production: The amplitude for the process 7 (g, X)j(q', A') —> tt + (p)tt~ (p'), 
is given by Al(A, A') = e 2 e fJ, (q, X)e"'(q', A')V)£, where the V^ v tensor can be decomposed into four 
Lorentz invariant tensor structures although by gauge invariance only two of them have non¬ 
vanishing contribution to the cross-section 

V£ v = A c {s , t, u)T lflv + B c (s, t, u)T 2[lv , 

Tiftu = - q^q'u ., T 2flv = 2sA l _ l A u - v 2 g^ u - 2u(q v A IJ - q'^ A„), (10) 

*To avoid confusion between the low-energy constants in xPT and S%PT the former are denoted by li while 
the latter by ii . The ii constants become li in the absence of S±. As is customary in xPT, the finite and scale 

independent terms of ii ( li ) are denoted by ii (li). 
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with s = (q + q ') 2 , t = (q — p) 2 , u = (q — p') 2 and := (p —The amplitudes A c (s,t,u ) and 
B c (s,t , u ) are analytic functions of the Mandelstand variables and are symmetric under crossing 
{t,u} -O- {u,t}- Comparison with the experimental data will be at the level of cross-section. The 
differential cross-section for unpolarized photons can be casted in terms of the helicity amplitudes 
corresponding to helicity changes A = 0, 2 respectively 


% = = \ H + + \ 2 + \ H +-\ 2 ■ ( n ) 

In terms of the amplitudes A c and B c they read 

H° + = A c + 2{4M$-s)B c , H%_ = 8(M ” ~ tu) B c . (12) 

At 0(p 2 ) there is no scalar contribution and the amplitude coincides with that of scalar elec¬ 
trodynamics. At 0(p 4 ) we have found remarkably many more diagrams than in yPT. Their 
evaluation is rather straightforward and the contributions can be conveniently cast in terms of 
two tensorial structures as 


_^( 4 ) _ e 2 A(s, ^ u)(se ■ e' — 2q ■ e' q' ■ e) + e 2 B(s , t, u ) 
which are related to those in (flUl) by 


, e- pe' ■ p' e • p't' ■ p 
e ■ e --— 


q-p 


q-p 


, (13) 


A c (s, t, u ) = 2A(s, t, u ) + 


B(s,t,u) ( 1 
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B c (s, t, u) = 


B(s,t,u) f 1 


4s 


+ 


Ml -t Ml — u 

1 


Ml - t Ml- 


(14) 


We have performed several checks on our full expressions: i) In the evaluation we have not fixed 
neither an specific gauge nor a system of reference and hence we are able to check explicitly gauge 
invariance in the results, ii ) All non-local divergences cancel when adding the full set of diagrams 
together with wave function renormalization, in) The polynomial divergences also cancel against 
the counter-terms determined in 0 2 75 = 76 = |(4 c 2 d — 1). iv) Once we shift the bare pion mass 
to the renormalized physical one the amplitude turns to be independent of Z\ and Z 2 . In view of 
these stringent checks, we trust our calculations of the matrix-elements. 

At this 0(p 4 ) order each of the above amplitudes can be split as 


A(s, t, u ) 


B(s , t, u ) 


-{2(2t 5 -l 6 )+G(s)} 


+ -p 2 A s (s, t, u) + 


CldCl-y 

F 2 




[2 ] + j±B a {s,t,u). 


(15) 


The explicit expressions for the A s (s,t,u), B s (s,t,u) and A 77 (s,t,rt) terms are gathered in the 
Appendix. The contributions in squared brackets correspond to %PT PI- Notice that corrections 
to B(s, t, u) at 0{p 4 ) are absent in yPT and only show up at higher orders JT4]. This confirms, as 
previously remarked in [7|, that the value for observables in SyPT at (D(p 4 ) lie within the 0(p 4 ) 
and 0(p 6 ) results in yPT. 


Results: To extract the value of the 77 Si coupling constant we have simultaneously fitted the 
experimental central values of the data for the processes n —> 717,77 —> 7r°7r° and 77 —> n + n~. 
For the latter we only take into account the data points in m near the two-pion production, 
yfs ~ 0.45 MeV, this removes to a large extent the KK effects. The data treatment of the 
former two experiments is described at lengthly in |7:. In all the procedure the only new free 
parameter, besides Ci 7 , at play with respect to those entering in [7] is the low-energy constant 
1 5 . We have generated a sufficient refined lattice for the set of constants, 5 x 10 6 points, in the 
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Figure 1: 77 —> tt + tt cross-section vs. the center-of-mass energy. Full line corresponds to the 
input values obtained in our fit as described in the main text msD, dashed line corresponds to yPT 
at 0(p 4 ). Dotted line corresponds to the results for the singlet obtained in JT6] where in addition 
to those we fitted the value of 24 — 4 = 0.0030. We used data values below 0.5 MeV for our fit, 
while the rest of the curve is just extrapolated. 


hyperplane defined by {c\d 1 M a , T', f/, 4, ci 7 jl3 with a priori flat distribution and computed their 
corresponding y 2 augmented function. Notice that we have treated all the coupling constants 
entering in the processes at the same footing, i.e. without imposing a priori any hierarchy, and 
nevertheless the output is consistent with the assumed counting power, |ci 7 | -C |ci/. The main 
result of this fit is given by 

cid = 0.26 +°;“ , l A = 2 4 - 4 = 0.0026±g;ggJt, 4 = 19.8^° 9 , 

M a = 553±?u MeV, V = 295/^ MeV, ci 7 = -0.012+°;°D® , (16) 


and is depicted in fig. |T] as the full curve for the 77 —> n + n cross-section. Fits to 77 —> 7r 0 7r° 
and 7r —> 7r7, not shown here, are similar to those obtained in (7|. The total y 2 , 0 j for the joint 
fit of all the three processes is 1 99 ' 4 . For comparison, the same fit but using yPT at 0(p 4 ) gives 
Xd o f = 3 65 4 ■ Notice that the finding concerning Ci 7 matches the short distance arguments that 
suggest a small two photon coupling 113. Errors in © correspond to the la deviations. It is 
worth emphasizing that the narrow thickness of the band in fig. [T| suggests that this experiment 
is not suitable to pin down the scalar mass and/or width. This statement is more evident if we 
compare our outputs for the singlet mass and width, (fl6l) , with those obtained in [16] , M a = 420 
MeV and P' = 286 MeV The latter, depicted as the dotted line in fig. [0 lies within the lcr 
deviation from the central value of (1161) . It is also worth emphasizing that a tiny variation in 
the fit, for instance including or not the data point at i/s = 0.395 MeV, changes the preferable 
set point that minimizes the data. 

The value of the combination of low-energy constants must be compared to those standard 
estimates obtained in [ 14 ] 2 1 5 — l 6 = 0.0028 and [l8j 2 1 5 — l e = 0.0031. Or to that extracted 
independently from the 7r + —> e + p e ') decay via the axial-vector-to-vector form factor ratio -p 1 


tNotice that is finite and scale independent. It can be expressed in terms of 7 quantities a.s f a = — „ l a = 

-9^(4-4). , 2 

tWe have rewritten the outputs of [I6l . where mass and width were defined as s a = (M a — ? '~rj > i n our 
convention s„ = M'i. — ir'M a . 
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2Z 5 — le = 0.0031 [T5]. The difference between those results and the corresponding one in (1161) 
gives an understanding of the effect of the singlet field in this combination of low-energy constants. 
As was expected from the beginning the contribution of the scalar singlet is mild in this process 
because it is mainly saturated by Vectors and Axials. 




Figure 2: Imaginary vs. real (Argand plot) parts of the Born subtracted 0(p 4 ) helicity amplitudes, 

_ Q _ (J 

H ++ {s,t = u ) and H . ( _ (s,t = u), as a function of the center-of-mass energy. The solid line is 

obtained using the central values in (1161) . The dotted curve is as above but setting to zero the 
electromagnetic coupling. Finally the dashed line corresponds to the yPT case. The dots signal 

_ (j 

the center-of-mass energy of the two-pion system in 100 MeV steps. Notice that H H does not 

receive any contribution from \PT at 0(p 4 ). Also the dashed line is indistinguishable from the 
solid line in this latter case. 


In fig. [2] we plot the imaginary vs. the real parts of the helicity amplitudes at t = u once the 
Born contribution is subtracted 

H++ := H° + - H° ++ , E%_ ,= H%_ - H c b+ _ . (17) 

It is evident that the electromagnetic correction is small and that at large energies there is a 
relatively large enhancement, with respect to the yPT, due to the inclusion of the scalar particle. 

As pass by we have also evaluated the dipole polarizabilities of the charged pion. This is 
obtained via the Compton scattering process y 7 r + —> 771 " + which is related to the pion-pair pro¬ 
duction by crossing symmetry s -o- t. Expanding m at the Compton threshold and using the 
input (1161) . we obtained 

( ai - /30.+ ^ (4.0±?;1, [6.0], {5.7}) x 10 " 4 fm 3 , 

(«! + 95 ( 0 . 012 ™ 4 , [ 0 ], {0.16}) x 10- 4 fm 3 , (18) 

where the numbers in square (curly) brackets stand for the standard yPT values at 0{p A ){0{p G )) 
respectively. As in the yPT case, it seems very hard to reconciliate the sharp discrepancy of (1181) 
with the most recent experimental result based on the radiative pion photo-production, yp —> 
7 K+n, ( ttl - PifJ = (11.6 ± 1.5 stat . ± 3.0 syst . ± 0.5 mod .) x 10 " 4 fm 3 ® . 
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Revisiting r(S'i — > 77 ): We are now in a position of finding the decay width of the scalar 
singlet to two photons. This was partially treated in [7] with the proviso that its direct coupling 
to photons was suppressed and the bulk of the contribution comes from the radiative process. 
Relaxing the above assumption and taking the 77 S 1 term into account we obtain 

nSi 77) = f^M a 3 |ci 7 - 8 c ld M ° M l M h (M 2 a )\ 2 - 0.1261q q 44 keV. (19) 

Notice that in the previous expression both terms, Born and radiative corrections, are of the 
same effective counting power. Analytically (1191) agrees with the Born approximation of [5] once 
we set ci d = 0. It is worth emphasizing the M 3 dependence in the above expression. This makes 
specially relevant the definition of the mass for a particle which width and mass are comparable. 
Had we used the convention in [16] our prediction for the central value of the scalar mass would 
have been approximately a 3% larger, or, equivalently, the sigma radiative width would increase 
a factor ~ 1.1. This can be accounted for as a source of systematic error. 

Owing to the smallness of (fTTJl) in comparison with the characteristic width of a conventional 
<50 resonance, for instance T(/ — ► 77) « 5 — 6 keV [21], we can conclude at the light of © that Si 
is mainly non-QQ. A comparison with other results that can be found in the literature is collected 
in table 1. One salient point is that ours is roughly a decade lower than the results obtained 
through dispersive calculations. 


Reference T(Si —»• 77 ) [keV] 


This work 


0 126 +u ' 34U 
U. 1 ZO — 0 044 


0 

m 

m 

m 

m 

m 

m 

m 

m 

[30] model A 

[501 model B 

m 


2.8 

1.68 ±0.15 
2.08 ± 0.20 

1.4 ~ 3.2 
3.08 
2.08 

1.7 ±0.4 
1.2 ±0.4 
4 

3.5 ±0.7 
2.4 ±0.5 

< 1 


Table 1: Comparison for T(Si —> 77 ) between different models. 

Aiming at a further theoretical interpretation we have checked whether this deviation w.r.t. 
the dispersive calculation can be assessed to a strong a —> KK coupling [25]. We have extended 
our analytical results to SU(3) and for a first and very crude estimation we took naively at face 
the values given in (1151) considering different ratios for the quantity 7^1 = By looking at 

the results, collected in table 2 , we may conservatively expect almost no sensitivity in the singlet 
decay width due to the presence of the strange quark mass. 

r a Kn 1 [universality] 0.8 [5] 0.37 [21] 0.62 [22] 

T(5i -> 77 ) [keV] 0T34 0T32 0T28 0T30 


Table 2: Comparison for an SU( 3) extension of T(S'i —> 77 ) using a naive extrapolation for the 
low-energy constants. 

In order to cross-check further our full approach we have computed the 77 — > 1 = 0 

s-wave phase shift in the threshold region. Those are related to the nn elastic scattering phase 

§ gairir stands for the obvious generalization of c] 4 . 
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shifts through Watson’s theorem. We have proceeded reconstructing the partial waves amplitudes, 
T[ , from the neutral and charged 77 —> tttt processes and through these the phase shifts. As is 
customary we express the 7r — 7 r elastic scattering result as the expansion in energy [52] 

5{ = arctan (Rel/) + 0{E & ) = <5, 7(2) + 5{ {A) + 0(E 6 ), 16 > s > 4. (20) 

At any time we bear in mind that the truncated chiral expansion becomes unreliable above y/s ~ 
450 MeV. In fig. [3] we have depicted our results for the central values (fl 6 l) adding for comparison 
the corresponding xPT ones. As is evident from the figure we obtain a remarkable improvement 
w.r.t. the yPT prediction and the agreement with the most recent data is rather good specially for 
the energy range 0.5 GeV < y/s < 0.7 GeV. We stress that there is no fit to these data and is just a 
prediction or a consistency check. This together with the fact that we reproduce the experimental 
data for the 7 r — 7 r scattering lengths, pion polarizabilities and the pion radii [7, let us to think that 
we have obtained a fairly good parameterization of the low-energy region containing the effects of 
the singlet state. 



y[7 (GeV) 

Figure 3: 77 —> 7 r + 7 r _ , 7 = 0 s-wave phase shift i n the low-energy region. The keys to 
experimental data are as follows: ■ = [33] , ♦ = [34|, A = [35], • = [36]. The blue dashed and red 
full lines correspond to the 0{p 4 ) result for yPT and SyPT respectively. The green dotted line 
denotes the nn elastic scattering Sq < ~ 2> . All curves agree at threshold, obeying Watson’s theorem. 


Once we have settled the consistency of the approach let us come back to the discussion on 
the quark content. The gluonium or four-quark scenarios are the most controversial scenarios to 
disentangle from our analysis. In order to do so we look at the relative weight between both terms 
in Considering just the direct coupling term one obtains the results in table 3. 


This work [3] [24] [25] [37] [38] [39] [40] Fit I [40] Fit II 

17,-.,, [keV] 0115^0? D.13 ±0.03 0.3 l).i(> 0.005 (>.03 -hi ] (jlj 0.024 ± 0.02.3 0.38 ± 0.09 


Table 3: Comparison of different results for the direct contribution to the decay width. 

Considering instead the decay via the rescattering process we get the results in table 4. The 

This work m mm 

[keV] 0T94^f| 2.7 ± 0.4 L89 2 


Table 4: Comparison of different results for the rescattering contribution to the decay width. 

difference with the 0.11 keV value found in [7] is due to the slightly bigger value of cu- Thus the 
initial mismatch with the dispersive calculation can be traced back to the rescattering term. In 
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particular the difference can have a twofold origin, see eq. m- o The constant C\d■ Although 
this constant is estimated at tree level its value is essentially upper bounded by the pion vector 
form-factor and the 1 = 0 s-wave phase shift <$g below the KK threshold. One expects that the 

( jyj- 2 \ _ 

-=r) ■ H) The G{M 2 ) function, 

or more generically pion rescattering effects. Notice that due to counting power contributions to 
Tg 1 _ >77 start already at 0(p 4 ) thus one expects higher order corrections of ~ 20% ~ 30%. To 
estimate these we have partially resummed a subset of higher order diagrams obtaining an increase 
of ~ 10% w.r.t. the central value in (fl9l) . Thus, although this result is incomplete seems to indicate 
that the numerical differences w.r.t. the dispersive results are hard to be asset to higher order 
corrections. 

On the other side one has to bear in mind that dispersive calculations are not free of uncer¬ 
tainties. Just to mention a few instances: i) The results seem to be very sensible to the matrix 
element parameterization above the KK threshold, ii) Only the s-wave component is kept at 
low-energy, iii) The result seems to be very sensitive to the actual value of the analogous of cid, 
i.G. Qct-k'k^ 0 For instance the differences between the value [24] [26j, given in table 1 and those 
found in i42j , 0.2 ~ 0.3 keV, are just due to the value of this coupling, iv) The approach, by 
analytical continuation, evaluates the matrix element deep in the complex plane. One has to keep 
in mind that the original embedding |43j is valid for point like particles which matrix elements are 
evaluated near the real axis. 

Comparing the central value for the direct and rescattering decay widths we learn that the 
relative weight between both terms in Cl is approximately 1 : 2 and that their interference is 
partially destructive. The relative smallness of the direct coupling in front of the radiative term, 
mediated via pion loops, can be interpreted as an indication of a dominant Q 2 Q 2 component in 
the nature of the scalar singlet. However this conclusion has to be taken cautiously as we have 
checked that for an increasing singlet mass the scenario can be reversed. Obviously all the above 
reflections are in the absence of mixing which can obscure this simple picture. In fact, this is 
neither strange nor new as similar conclusions are supported by QCD sum-rules [2j, lattice QCD 
calculations [33] and large N c scaling arguments [3Sj . The novelty of our approach resides in that 
this finding is encoded in the low-energy regime and an effective approach suffices to capture it. 

Conclusions: We have found an estimate to the scalar to two photons decay width using low- 
energy data. The fact that the preferred point is attained for |ci 7 | < 1 but not vanishing signals 
the presence of the canonical anomaly irn- Making use of the central values of © together 
with m, and R = | for consistency, we obtain Fs 1 ~ —2.3 F n , to be compared with the QQ 
and the pure glueball results: Fs 1 = —F^ and Fs 1 ~ —5 F„ respectively 0. This fact reinforces 
our conclusions about the tetraquark nature of the scalar meson as derived from its coupling to 
two-photons ©• 

Adopting the most optimistic attitude, taking into account higher order resummations, SU (3) 
extensions and the sensitivity of the results on c\d and M a the central value in (fH)l) can be pushed 
up to 

r(Si 77 ) « (0.3 - 0.4) keV. (21) 

This agrees with other effective approaches, studies of the low-energy data using a Breit-Wigner 
cross-section or studies of the 77 — > ir°Tr° cross-section assuming a scalar dominance I Hi . There 
is however a mismatch of a factor 4 ~ 5 w.r.t. the dispersive approaches. 

Concerning the possible sources of difference w.r.t. the dispersive calculations two comments 
are in order: i) We have an analytic expression for the rescattering piece at 0{p 4 ). We remind that 
within the effective framework unitarity is only satisfied perturbatively contrary to the dispersive 
approach where unitary is enforced by construction and use of high-energy data is taken into 
account. This is also the main reason underlaying the small deviation from the scattering phase 
shifts data below 0.5 GeV as in the standard case [35] . ii) Concerning the second source, the 
coupling cid , it has a more controversial status. Being a tree level constant we have found it 


* Not to be confused with the energy independent generalization of c ,,/ used in Table 2. 



essentially through processes in which its role enters at the radiative level. Due to the sensitivity 
in the dispersive approach to the value of g a mr(s) it would be interesting to have a constraint on 
cid in processes where it plays a dominant role even at leading order. 

We stress that only low-energy data were used in our approach. 
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Appendix 

We have gathered in this appendix all the relevant information concerning the rational functions 
and integrals appearing in the amplitudes (TTSl) together with the diagrams that describe the process 
77 —*■ 7 T+ 7 T - , see fig. ( 0 ]). In the calculation we used dimensional regularization in the MS scheme. 
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(9) 
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) 


(12) 

(13) 



(23) 




Figure 4: Feyman diagrams for the process 77 n + n~. The keys to the states are as follows: 
wavy lines=photons, full lines = pions and double full lines = scalar singlet. Diagrams from (1)- 
(13) denote the s-channel and they contribute to the e • e' component of the amplitude. The t 
(u)-channel, diagrams (14)-(25), contribute to the transverse components of the amplitude. 


The short hand notation of the amplitude can be casted in terms of the finite part of the one-, 
two-, three- and four-point scalar functions as 

^ 0 ( 5 ,^,!^) — [4 7 r [4(7 , /Ci(s,t,u) — J tt a () , A 2 (s, t, u) — J 77 a (t) , 

/C 3 (s, t, u) = G{s ), /C 4 (s, t, u ) = Co (a, M=, M 2 a ) , 

fCs(s, t, u ) = Co (0, t, Ml Ml Ml Ml) , 

/C 6 (s, t, u ) = Do (0,0, Ml Ml s , t, Ml Ml Ml M 2 ) , (22) 


where the Co and Dq functions are the ones introduced in 1491 and the overline indicates that they 
incorporate the yJ-j factors, as the J and G functions do. In particular 


G{s) 


1 

167T 2 


2M 2 C o (0,0,s,MlMlM 2 ) 


(23) 
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We can split the amplitudes in terms of a polynomial piece and a dispersive one as 


A s (s,t,u) = 
B s (s,t,u) = 


P A (s,t,u) + ^2 U A i{s,t,u)ICi(s,t,u) 
?—o 
6 

P B (s,t,u) + U B i(s,t,u)ICi(s,t,u) 


2—0 


+ [t 0 u] , 
T [t -O- u] , 


(24) 


where the Kf s correspond to the scalar loop functions and U's are rational functions of the masses, 
scalar width and Mandelstand variables, with v = t — u. The terms contributing to the A s (s, t, u) 
contribution are 


P A (s,t,u) = 


s (Ml - 2 Ml) 2 (2 Ml - AMI + a) 


107T 2 (Ml — tu) 2 

U A0 (s,t,u) = 0, 

4 (Ml - 2 Ml) 2 (M^s + 2Mltu - 2AT® + stu ) 


U A i(s,t,u) = 


U A2 (s,t,u) = - 
U A3 (s,t,u) = 


s (t — Ml) (Ml - u) (tu - Ml) 
8 1 (Ml - 2 Ml) 2 (Ml - u) 
s(t-Ml)(M$-tu) ’ 

1 


U A4 (s,t,u) = - 
U A5 (s,t,u) = - 
U A6 (s,t,u) = - 


8s (Ml - tu ) 2 L 16 ^ ^ 10sM ' + fi2 + ^ 2M2 ( fi2 - v2 ) +s ( fi2 - ^ 
+ 8s 3 Ml + 16s 2 Ml - 6AMls 2 - AMl (2As 2 Ml + AsM 2 a (s 2 + v 2 ) + s 4 - t/ 4 )] 

-4(s-2M 2 ) 2 ^, 

s 

2Ml (Ml - 2Ml) 2 (—2 Ml + 1 2 + u 2 ) 


(Ml — tu) 2 

8 M 2 (M 2 - 2M 2 ) 2 (M 2 - t) (-Ml + Af 2 + 1) 
(Ml — tu) 2 

AMl (Ml - 2M 2 ) 2 


s 2 M 4 + Ml (s 2 - 12t 2 ) - 2s 2 Ml (Ml + t) 


s (Ml — tu) 2 

+2M 2 f 2 (3s + 4i) - 2Ml(s - 4t) - 2M® - f 2 (s 2 + 4sf + 2t 2 )] , 
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and the terms contributing to the B s (s , i, u) amplitude read 


PB(s,t,u) - 

U B0 {s,t,u) 

U B i(s,t,u) 

U B2 (s,t,u) 
U B z{s,t,u) 
U B4: {s,t,U) 
U B5 {s,t,u) 
U B6 (s,t,u) 


- (M| - 2 Ml) 2 

16tt 2 M2 (4M 2 - M2) (M4 - tu) 2 


[4M 2 (M 2 s 3 + Af4 (_ 6 s 2 + 1 2 + u 2 ) + 8 M°s 


— 5M® + tu ( t 2 + tu + u 2 )) + 2sM% (i — M 2 ) (M 2 — u) 
+ s (s — 12M 2 ) M 2 (t — M 2 ) (Ml-u)] , 

8 M|(M 2 - 2 M 2 ) 2 
: -5M 2 M 2 + M 4 + 4M 4 ’ 


4 (M 2 - 2M1) 2 r s , , 

Ml {AMI - M 2 ) (Ml - t) {Ml - u ) {Ml - tu) L w ^ a ’ 
+Mltu (2 sM'l - 3 tu) - 7Mlstu + fVM 2 - 2M%tu + 5M 30 ] , 
8 (M 2 - 2 Ml) H (u - Ml) 

(t ~ Ml) {tu — Ml) ’ 

s (Ml - 2Ml) 2 (t - Ml) (Ml - u) (2Ml - 4M 2 + s) 

Ml {Ml - tu ) 2 

2 (Ml - 2 Ml) 2 (t - Ml) (u - Ml) (-2MI + t 2 + u 2 ) 
(Ml-tu) 2 ’ 

8 (Ml - 2 Ml) 2 (t - Ml) 2 (Ml - u) 

(Ml-tu) 2 

4 (Ml - 2Mp 2 (t - Ml) (Ml - u) 

{Ml — tu ) 2 

+ Ml{s + At) - 2Ml\ . 


(_ Ml+Ml+t) . 

[- 2Ml (sMl +t 2 )+s(t- Ml) 


Finally, the amplitude proportional to the direct <777 coupling ci 7 (fT5l) is 


I 77 (s, t, u) = (| - M 2 ) F(s). 


(25) 


Notice that in the above expressions we have used a Breit-Wigner representation to regularize the 
propagator of the scalar particle 


F(s) 


1 

s- Ml + iM a r ' 


(26) 


This can be, at first sight, slightly controversial. The main two arguments to use this parametriza- 
tion are: i) as in all the processes studied in J7] in this work the propagator enters in the highest 
radiative order, thus differences between parametrizations would be reflected at least at 0{p 6 ), 
beyond our scope. This would drastically change in the case of studying it — -k scattering where 
already the scalar propagator enters at lowest order, ii) In this line, we have recovered the results 
in |i7], within the la band, using a different parameterization 1501 : 


F(s) 


1 

s - Ml + iM a r{s) ’ 


with T(s) 


s - s 0 
Ml - so 


3/2 


r 0 . 


(27) 
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